Color confinement and dual superconductivity of the vacuum. III. 



O 

o 

(N 



> 

in 
o 
o 
m 
o 



I 

Oh, 



X 



Dipartimento 



J.M. Carmona^'i, M. D'Elia'''^, A. Di Giacomo'=''''^ B. Lucini'^^^^ q Paffuti^''^-^ 
" Departamento de Fisica Teorica, Universidad de Zaragoza, 50007 Zaragoza, Spain 
Fisica dell'Universita di Genova and INFN, Sezione di Geneva, Via Dodecaneso 33, FI6I46 Geneva, 

Dipartimento di Fisica dell'Universita, Via Buonarroti 2 Ed. B, 1-56127 Pisa, Italy 
INFN sezione di Pisa, Via Vecchia Livornese 1291, 1-56010 S. Piero a Grado (Pi), Italy 
Theoretical Physics, University of Oxford, 1 Keble Road, 0X1 3NP Oxford, UK 
^ e-mail address: jcarmona@posta.unizar.es 

^ e-mail address: delia@ge.infn.it 
^ e-mail address: digiaco, paffuti@df.unipi.it 
* e-mail address: lucini@thphys.ox.ac.uk 

It is demonstrated that monopole condensation in the confined phase of SU{2) and SU{3) gauge 
theories is independent of the specific Abelian projection used to define the monopoles. Hence the 
dual excitations which condense in the vacuum to produce confinement must have magnetic U{1) 
charge in all the Abelian projections. Some physical implications of this result are discussed. 

PACS numbers: 11.15.Ha, 12.38. Aw, 14.80.Hv, 64.60. Cn 



Italy 



I. INTRODUCTION 

This paper is the third of a scries, in which we study 
the dual superconductivity of the QCD vacuum as a 
mechanism for confinement of color. In the first two pa- 
pers dHI] we have detected condensation of monopoles in 
the confining phase by means of a disorder parameter (fi) , 
which is the vacuum expectation value of a magnetically 
charged operator, /i. (/i) was determined by simulating 
the theory on a lattice: it is non zero in the confining 
phase, and tends to zero at the deconfining transition, 
above which it vanishes. 

The connection of (/i) to confinement was proved by 
a quantitative determination of the critical indices and 
of the critical coupling. In H), SU{2) gauge theory was 
studied, while Q was devoted to SU{3) with similar re- 
sults. 

Magnetic charges in gauge theories are defined by a 
procedure known as Abelian projection j^]: to every local 
field $ belonging to the adjoint representation {Nc — 1) 
f/(l) fields can be associated, and to each of them a con- 
served magnetic charge. In fact there exists a functional 
infinity of monopole species, A^c — 1 for each field <I>, which 
in principle can condense in the vacuum and confine the 
corresponding U{1) electric charge by the dual Meissner 
effect. It is not known a priori if monopole condensations 
in different Abelian projections are independent phenom- 
ena. 

The indication obtained in |^,|| by analysis of a number 
of different choices of <f> was that all of them show the 
same behaviour, so that they are equivalent to each other. 

The possibility that in some way all the Abelian projec- 
tions could be physically equivalent was first advocated 
m Ref. §. In this paper we add strong evidence for that 
equivalence. 

In order to explain what we do, let us first recall how 
magnetic charges are associated to any field $ in the ad- 
joint representation. We shall do it for SU (2) to simplify 
notations; extension to SU{N) only adds formal compli- 



cations. 

Let $(x) be a field in the adjoint representation (color 
vector), and let <I>(x) be its color orientation, 



<l(x) 



$(a;) 



(1) 



is well defined, except at zeros of $(a;). 
Define a gauge invariant field strength Ffj_^{x) 



(2) 



where G^^, = d^,A^ - d^A^ + gA^ A A^, is the gauge 
field strength and Z3^$ = (9^ + .gyl^A)$ is the covariant 
derivative of $. 

Both terms in the right-hand side of eq. (||) are sepa- 
rately gauge invariant and color singlets: their combina- 
tion is chosen in such a way that bilinear terms A^^A^, 
y4^$ and 71^$ cancel. Actually, by simple algebra, 

F^, = $ • {d^A, - d,A^) - i A dA) • $ . (3) 

If we transform to a gauge in which = constant 
in space-time, the last term cancels and 



(4) 



is an Abelian field strength. Such a gauge transforma- 
tion is called an Abelian projection. It is in general a 
singular transformation which exposes monopoles at the 
sites where ${x) — 0. 

If = ^€fj,,ypaFP°' is the dual field to F^i^, one can 
define a magnetic current 



d^F* ■ 



(5) 



ifj, is zero if Bianchi identities hold, but can be non zero 
in compact formulations in terms of parallel transport. 
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like lattice formulation ||7|. In any case it follows from 
the antisymmetry of F*jj that 

d^], - . (6) 

In the dual superconductor view of color confinement 
the symmetry (|^) is expected to be realized a la Wigner in 
the deconfined phase, and to be broken d la Higgs in the 
confined phase. An operator fi which carries magnetic 
charge can provide a disorder parameter to discriminate 
between the two possibilities. Such an operator was de- 
veloped and tested in Refs (|-|l^,§J|| . 

What was found in was that there is indeed dual 
superconductivity in a number of Abelian projections. 
As explained in detail in Sect. II, the full identification of 
the projected gauge requires to go to the gauge in which 
${x) - A is diagonal in color indices (A are the generators 
in the fundamental representation), with a fixed order of 
the eigenvalues. 

One can diagonalize $(a;) - A up to the ordering of the 
eigenvalues, choosing it randomly, and still define an op- 
erator fi which creates a magnetic charge in that Abelian 
projection. We show in this paper that the correspond- 
ing disorder parameter behaves exactly in the same way 
as the one with ordered eigenvalues. 

We can even define a completely random Abelian pro- 
jection, in which we do not diagonalize any operator $, 
but we take, e.g. for SU{2), $ — 0-3, the nominal 3 axis 
used in the simulation, and define the corresponding fi. 
Again we find that fi defined in this way behaves exactly 
in the same way as those defined in Refs. ^,|| and scales 
with the same critical indices. 

The above Abelian projections are kind of an average 
over a continuous infinity of Abelian projections, and the 
result demonstrates, beyond any doubt, the complete in- 
dependence of dual superconductivity from the choice of 
the Abelian projection. 

Our results are compatible with Ref. where our 
disorder parameter in the random gauge was computed 
by Schroedinger functional techniques. 

There has been in the literature in the last years the 
idea that monopoles defined by a particular Abelian pro- 
jection (the maximal Abelian projection) are more rele- 
vant than others to confinement . We will discuss 
this issue in Sect. V, where we draw conclusions from our 
results. 

In Sect. II the construction of the disorder parame- 
ter (/i) will be recalled, to define the Abelian Projec- 
tion with Random Ordering (APRO) and the Random 
Abelian Projection (RAP). In Sect. Ill the numerical al- 
gorithms used will be discussed. The results will be de- 
scribed in Sect. IV. Sect. V will close the paper with 
conclusions. 



II. DISORDER PARAMETER 

In this Section we recall the definition of the disorder 
parameter for confinement. 

Let O be an operator which transforms in the adjoint 
representation of the gauge group, i.e. 

with A" the generators in the fundamental representa- 
tion. The Abelian projection technique ||] prescribes to 
fix the gauge by a gauge transformation in such a way 
that 

Ogf =G^OG = dmg{oi,...,0N) 

with oi < 02 < ... < On . (8) 

After Abelian projection, there is still a U{1)^^^ gauge 
freedom left, since a transformation of the form 

C! = diag(e'^"S...,e^-") , J]^. =0 (9) 

does not change the gauge fixing condition, Eq. (||). 

After Abelian projection, the gauge variables of 
SU{N) are divided in two sets: the photons (the N — 1 
neutral fields under the residual U{1)^~^) and the gluons 
(charged fields with respect to the residual symmetry). 
Abelian magnetic monopoles can arise at points where 
two eigenvalues of O are degenerate [|]. 

Condensation of Abelian monopoles defined by 
Abelian projection has been demonstrated numerically 
in Ref. g for SU{2) and in Ref. § for SU{3). This 
has been done by constructing an operator magnetically 
charged in a given Abelian projection and by studying the 
behaviour of the vacuum expectation value of that oper- 
ator across the phase transition at finite temperature. In 
the language of Statistical Mechanics we call that opera- 
tor a disorder operator and its vev a disorder parameter, 
the terminology being that the weak coupling (decon- 
fined) phase is the ordered phase. The construction can 
be done in different Abelian projections: in Refs. a 
number of Abelian projections have been studied, and for 
all of them it was found that indeed monopoles condense 
at low temperature, while the corresponding magnetic 
symmetry is implemented a la Wigner at high tempera- 
ture. Moreover, the disorder parameter scales with the 
correct critical indices in the critical region and is inde- 
pendent of the choice of the Abelian projection. These 
results suggest that the observed behaviour of the dis- 
order parameter is generally independent of the Abelian 
projection and of the Abelian operator chosen. 

Let us review the construction of the disorder param- 
eter. We introduce a time-independent external field 

$,(n,y) =Ge'^'''("~^'S')Gt , (10) 

where G is the gauge transformation that diagonalizes 
the operator O according to Eq. (^, & is the discretised 
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transverse field {i.e. V-6 = on the continuum) generated 
at lattice spatial point n by a magnetic monopole sitting 
at y and T is a generator of the Cartan subalgebra. 

Let Ufii, be the Wilson plaquette defined in the usual 
notations as 

U^,,ix) = U^ix)U,{x + fl){U,,ix + P)y{U,ix))^ , (11) 

where x = {n,t). We introduce a shift Uio{n,0) 
Uio{n,0) by inserting the external field $i(n + i,y) in 
the path ordered product at time zero, as follows: 

U^o{n,0) = U^{n,0)^,{n + i,y)Uo{n + i,0) (12) 

and we define U^v{x) = Ufj,u{x) elsewhere. 

The Wilson action for SU (N) gauge theory is 

where the sum extends over all the lattice points and 
directions. 

By replacing in the previous equation the standard pla- 
quette with the modified plaquette U^y{i) we obtain the 
"monopole" action 

5M(y,0) = (^1 - ^ (ii^^ix) + (?7^.(x))t)) . (14) 

is given by 

0)) = ' 'V,':,, , (15) 



2. We do not perform the Abelian projection, i.e. we 
take 



The disorder parameter introduced in 

/(PC/)e-^»-^(S'o,o) 



J{VU)e-S ' 

where the functional integral of e~^ is taken with peri- 
odic boundary conditions and the integral of e^^''' with 
C*-pcriodic boundary conditions ||lj,|l^ , 



U,(n,t = Nt)^U*[n,t^Q) , 



(16) 



Nt being the temporal extension of the lattice and U* 
being the complex conjugated of Ui. 

To study more in detail the dependence on the pro- 
jecting operator, we will modify the definition of $ as 
follows: 

1. We choose a projecting operator and we diagonalise 
it but without fixing the order of the eigenvalues, or 
better with the order of the eigenvalues randomly 
chosen: 



(17) 



where P is a random N x N permutation matrix. 
This corresponds to a sort of average of fx over the 
class of operators differing from O on each point 
by the order of the eigenvalues. We refer to this 
case as Abelian Projection with Random Ordering 
(APRO). 



$,(n,y) ^e'^'''^"-^'^^) 



(18) 



This is equivalent to a sort of average of over all 
the possible Abelian projections. We refer to this 
case as Random Abelian Projection (RAP). 

As discussed in Ref. a direct computation of (^} 

with Monte Carlo techniques is problematic, because this 
quantity has large fluctuations, being the exponential of a 
sum over the physical volume. A more convenient quan- 
tity to study in numerical simulations is 



d 

-gp log(Ai) = {S)s - (S'm)sa 



(19) 



p is the difference of two average actions, the 
Wilson action and the modified action Sm (the 
latter being averaged with the modified measure 
((X'C/)e-'5")/(/(X'[/)e-'5")). p has smaller fluctuations 
and contains all relevant information. The value of (p) 
is related to p by the relationship 



(M)=exp( 1^ p(Wj . (20) 



III. GAUGE FIXING AND SIMULATION 
ALGORITHMS 

We have determined the temperature dependence of p 
for SU{2) and S'f/(3) pure Yang-Mills theories, for both 
definitions (^^ and of $ on an asymmetric lattice 
X Nt with Nt<^Ns. 

For both definitions of $, the simulation of the Wilson 
term, (5)5, has been performed on a lattice with peri- 
odic b.c, by using a standard mixture of heatbath and 
overrelaxed algorithms. 

As for the APRO case, we have chosen the Polyakov 
line as the operator to identify the Abelian projection, 
following the definition in Eq. (31) of j|] and Eq. (19) of 

, with the only difference that at each spatial point the 
ordering of the eigenvalues is selected randomly among 
the possible different permutations Up (rip — 2 (6) for 
SU{2) {SU{3)) pure gauge theory). This effectively 

corresponds to averaging over Up " different definitions 
of the Abelian projection. The Abelian generator 
(F^ = A^/2, with A* the Gell-Mann matrices), has been 
chosen to define the monopole field for the SU{3) case. 
We use C* boundary conditions in time to compute 
(Sm)sm i'^ Eq. (^. In this case, as explained in Ref. 
[Q , it is not possible to use a standard heatbath or over- 
relaxed algorithm to simulate the modified action, since, 
6.17. in the case of Polyakov projection, the change of any 
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temporal link induces a non linear change in the modi- 
fied action. So we have performed simulations by using 
a mixed (heatbath + overrelaxed) algorithm for the up- 
date of spatial links, and a metropolis algorithm for the 
update of the temporal links. 

In the case of the Random Abelian Projection instead, 
as it appears in Eq. (JTsj), one does not need any gauge 
fixing to define the monopole field. As a consequence, 
the change of any link always induces a linear change 
in the modified action. Therefore we have used a stan- 
dard (heatbath + overrelaxed) algorithm in this case. In 
the SU{3) case, the Abelian generator has been used 
again to define the monopole field. 



IV. NUMERICAL RESULTS 

The phase transition is known to be second order for 
SU{2), weak first order for SU{3) As usual we shaU 
speak of critical indices in both cases, meaning for SU (3) 
effective critical indices at small values of (1 — T/Tc), but 
not too small. 

In Refs. jJJ^ it was shown that the critical indices 
of the confinement transitions for SU (2) and SU{3) did 
not depend on the particular type of Abelian projection 
used to define the monopole condensation. Here we will 
show numerically that the critical exponents, and also the 
value of the critical coupling, are the same even for the 
Random Abelian Projection and the Abelian Projection 
with Random Ordering, both for SU{2) and SU{3) 

The critical behaviour of the disorder parameter (/i) 
is governed by an exponent 6. For finite lattice sizes 
{N^ X Nt), finite size scaling states that 



(21) 



where a and ^ are respectively the lattice spacing and the 
correlation length of the system. 
Near the critical point, for /3 < f3c, 



ecx(/5c-/3)- 



(22) 



where v is the corresponding critical exponent. In the 
limit Ns ^ Nt and for a/^ ^ 1, i.e. sufficiently close to 
the critical point, 



or equivalently 
P 



N. 



'nI'" {(3c - P) 



(23) 



(24) 



The ratio p/Nl^'^ is a universal function of the scaling 
variable 



We will use the known values of fic and v of SU{2) and 
SU{3) pure gauge theories to see that scaling holds with 
the present data. In order to obtain the critical exponent 
S, we use an expression equivalent to Eq. (|2^), 



From here we get 



(/3, - PYt{x) 



N. 



l/z. 



(26) 



(27) 



To obtain S we need additional assumptions on the un- 
known scaling function J-{x) . We will see that fits of good 
quality are obtained with the simple parametrization 



P 



= c, 

X 



(28) 



where C is a constant term. This form is suggested by 
the fact that when 2: ^ 0, both T{x) and its derivative 
should go to a constant. 



A. S'C/(2) gauge theory 

1. The Random Abelian Projection 

The quality of the scaling, Eq. (|2|), for SU{2) in the 
RAP case can be seen in Fig. |^. Here we used the known 
values of /?c = 2.2986 and 1/ = 0.63 0. The curve in 
the figure corresponds to the best fit to Eq. (|28|). We get 
S — 0.19(5), with a x^/d.o.f. ~ 1.5, in good agreement 
with the value obtained in the plaquette and Polyakov 
gauges in Ref. g, (5 = 0.20(8). 




N,""(|3c-P) 



FIG. 1. Quality of scaling in the RAP case for 5';7(2). 



= iVy^' (/3c - /3) . 



(25) 
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and 1/ = 1/3 psf . The curve in the figure corresponds 
to a best fit to Eq. (|2^), modified by including the 
term ^'(A^s), which gives the value 6 = 0.50(3), with a 
xVd.oi. = 3.2. 

The quality of the fit improves considerably if one in- 
cludes a second term in the expression for ^'(TVs): 



(31) 



Then the value of 6 remains the same, while x^/d.o.f. ~ 
0.7. This demonstrates the importance of the finite size 
efi^ects in SU (3) gauge theory. The fit is shown in Fig. |^. 



1 2 3 4 5 

N."^'(|3c-P) 

FIG. 2. Quality of scaling in the AFRO case for SU{2). 

2. The Ahelian Projection with Random Ordering 

Results obtained in the AFRO case are shown in Fig. 
where again known values of j3c and v have been used. 
The curve in the figure corresponds to the best fit to 
Eq. (|8|), which gives 5 = 0.20(6), with a xVd-o.f. 1.1. 
The agreement with the results obtained in the RAF case 
and in the plaquette and Folyakov gauges in Ref. [Q is 
very good. 

B. SU{'i) gauge theory 

The confinement transition in pure SU (3) gauge the- 
ory is a first order transition ||l^]. One therefore expects 
a pseudocritical behaviour, with v = 1/3, that is, the 
inverse of the number of spatial dimensions. As it was 
remarked in Ref. ||^, the scaling relation, Eq. (p4|), has to 
be modified in this case to include finite size violations 
to scahng. 



N 



^ = / [N]!'^ {(3c - P)\ + ^{Ns) 



(29) 



where '^{Ng) parametrizes these effects. A simple as- 
sumption is 



valid up to 0{l/N'^). 



a 



(30) 



1. The Random Ahelian Projection 

Fig. ^ shows the scaling behaviour expressed by 
Eq. (|29|), where we have taken for '^{Ng) the form (30). 
We used as an input the values Pci^t = 4) = 5.6925 



i -0.05 - 




5=0.50(3) 
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FIG. 3. Quality of scaling in the RAP case for SU{3). 




6=0.50(3) 
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N=16 
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N=32 





100 200 300 

N,""(Pc-P) 



FIG. 4. The same as in Fig. |, using Eq. (|l|) for *(iV., 
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2. The Abelian Projection with Random Ordering 

In the APRO case we have obtained a good scahng 
behaviour as well, as shown in Fig. ||. 

However in this case a fit to Eq. ( p9| ) with the func- 
tion ( |30| ) has a very bad x^/d.o.f. (of order 16) . 

Also in this case the use of the expression is essen- 
tial: the best fit, shown in the figure, gives 6 = 0.43(3) 
with x^/d.o.f. ^ 1.4. The value of S is nearly compatible 
with the one obtained in the RAP case, S ~ 0.50(3). The 
result obtained in Ref. @ is (5 = 0.54(4). 



one single Abelian projection does not confine the U{1) 
neutral particles belonging to the adjoint representation, 
which would instead be confined in other Abelian projec- 
tions. 

Analysis of the Zn vortices could give some hints and 
investigation has been started in that direction |2^. We 
think that the problem is still open. 
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FIG. 5. Quality of scaling in the APRO case for SU{3), 
using Eq. {B^ for ■^{Ns). 



V. CONCLUSIONS 

We have produced further and compelling evidence 
that monopole condensation is independent of the 
Abelian projection used to define the monopoles. 

If the idea of duality is correct, the non-local excita- 
tions which are expected to be the fields of the dual de- 
scription of QCD and weakly interacting in the confined 
phase, should have non zero magnetic charge in all the 
Abelian projections. This is a very important symmetry 
property, which can help in identifying them. 

There has been a number of papers in the literature 
of the past years, claiming that the fundamental fields 
of the dual description are the monopoles defined by the 
maximal Abelian projection. The claim that monopoles 
defined by the maximal Abelian projection could be 
the dual excitations does not look in good shape af- 
ter the quantitative attempts to construct the dual the- 
ory, which go beyond the initial empirical observation of 
Abelian dominance [p^ . If it were true maximal Abelian 
monopoles should be magnetically charged in all Abelian 
projections. This does not seem to be plausible, since 
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